Accepted for publication in the Astrophysical Journal Letters 
Preprint typeset using LfTgX style emulateapj v. 04/03/99 



CONSTRAINTS ON REIONIZATION FROM THE THERMAL HISTORY 
OF THE INTERGALACTIC MEDIUM 1 ' 2 

Tom Theuns 3 , Joop Schaye 4 , Saleem Zaroubi 5 , Tae-Sun Kim 6 , Panayiotis Tzanavaris 3 

and Bob Carswell 3 
Accepted for publication in the Astrophysical Journal Letters 

ABSTRACT 

The temperature of the diffuse, photo-heated intergalactic medium (IGM) depends on its reionization history 
because the thermal time scales are long. The widths of the hydrogen Lya absorption lines seen in the spectra 
of distant quasars that arise in the IGM can be used to determine its temperature. We use a wavelet analysis of 
the Lya forest region of quasar spectra to demonstrate that there is a relatively sudden increase in the line widths 
between redshifts z«3.5 and 3.0, which we associate with entropy injection resulting from the reionization of 
He II. The subsequent fall-off in temperature after z ~ 3.5, is consistent with a thermal evolution dominated by 
adiabatic expansion. If, as expected, the temperature also drops rapidly after hydrogen reionization, then the high 
temperatures inferred from the line widths before He II reionization imply that hydrogen reionization occurred 
below redshift z = 9. 

Subject headings: cosmology: observations — cosmology: theory — galaxies: formation — intergalactic 
medium — quasars: absorption lines 



1. INTRODUCTION 

Neutral hydrogen in the intergalactic medium (IGM) along 
the line of sight to quasars at redshifts z < 6 produces hundreds 
of Lya absorption lines. The fact that not all flux is absorbed 
(i.e., the absence of a 'Gunn-Peterson' trough, Gunn & Peter- 
son 1965) requires that the universe be ionized to a level far 
higher than can be attributed to residual ionization from recom- 
bination. At lower redshifts, z <; 3, observed stars and quasars 
produce enough ionizing photons to explain the high levels of 
ionization, but the nature of the sources responsible for convert- 
ing most of the IGM from neutral to ionized remain uncertain, 
as does the epoch of reionization (e.g., Barkana & Loeb 2001). 

The observed mean flux decrement Da blueward of the 
quasar's Lya emission line increases with the redshift of the 
quasar, both because the intensity of the ionizing background 
radiation decreases above z = 4 (e.g., McDonald & Miralda- 
Escude 2001) and because the mean density p of the universe - 
and hence the neutral fraction x oc pT~° J /Tm for fixed values 
of the photo-ionization rate Thi and temperature T - increases. 
Recently, Becker et al. (2001) and Djorgovski et al. (2001) ob- 
served a sudden increase in Da in the spectra of redshift z ~ 6 
quasars discovered by the SLOAN digital sky survey. Such a 
sharp rise has been predicted to mark the transition associated 
with a sudden epoch of reionization (e.g., Cen & Ostriker 1993; 
Gnedin 2000). Similarly, a sudden increase in the He II opacity 
has been detected around z ~ 3 (Reimers et al. 1997; Heap et al. 
2000; Kriss et al. 2001), associated with helium reionization. 

Another way to study the reionization history of the IGM is 
to investigate its thermal evolution. Because its cooling time 
is long, the low-density IGM retains some memory of when 
and how it was reionized (e.g., Miralda-Escude & Rees 1994; 



Hui & Gnedin 1997; Haehnelt & Steinmetz 1998). The com- 
bined effects of photo-ionization heating and adiabatic expan- 
sion introduce a tight temperature-density relation in the un- 
shocked IGM, which can be approximated by a power-law 
T = T()(p I py~ l for densities around the cosmic mean (Hui & 
Gnedin 1997). A change in these parameters influences the 
shapes of the Lya lines, because thermal broadening and Jeans 
smoothing determine the line widths (e.g., Theuns, Schaye & 
Haehnelt 2000). Schaye et al. (1999) used hydrodynamical 
simulations to demonstrate that one can accurately calibrate the 
relation between the minimum line width (b) as a function of 
column density (A/hi) on the one hand, and the underlying T-p 
relation on the other. Schaye et al. (2000) applied the method 
to observations in the redshift range 2.0-4.5 and found that Tq 
peaks at z ~ 3, which they interpreted as evidence for the reion- 
ization of He II. Ricotti, Gnedin & Shull (2000) used pseudo 
hydrodynamical simulations and found a similar temperature 
increase, albeit only at the 0.5er level. McDonald et al. (2001) 
found no evidence for temperature evolution, but their analy- 
sis neglected the important temperature dependence of Jeans 
smoothing. Finally, the analysis of Zaldarriaga et al. (2001) 
neglected hydrodynamical effects all together. 

Here, we provide new evidence for a relatively sudden in- 
crease in 7b between redshifts z ~ 3.5 and 3.0, using a new 
method based on a wavelet decomposition of the absorption 
spectrum. We then use the measured values of 7b at higher 
redshift to constrain the epoch of hydrogen reionization zh 
and find that the data require zh < 9 for any reasonable value 
(To <J 6 x 10 4 K) of the hydrogen reionization temperature. 
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2. HELIUM REIONIZATION 
2.1. Wavefef analysis 

Our analysis uses wavelets to characterize line-widths. By 
picking an appropriate wavelet scale (we used ~ 15 km s" 1 ), 
we find that the amplitude A of the wavelet anti-correlates with 
the widths of the lines - and hence the temperature 7b of the ab- 
sorbing gas, (A) oc r _1 . By examining the statistics of A along 
the spectrum, we can look for changes in temperature in an ob- 
jective way, given that the wavelet decomposition is unique. 
Full details can be found in Theuns & Zaroubi (2000) and The- 
uns et. al (2002), here we give only a brief summary of the 
underlying reasoning. 

To investigate whether a region of given size V of the QSO 
spectrum has an unusual temperature, we compare the cumula- 
tive probability distribution Cy(A) of the wavelet amplitudes in 
that region, with C(A) for the spectrum as a whole. If the region 
is unusually hot, it will tend to have very few large wavelet am- 
plitudes, and hence the maximum difference A(V) = Cy-C will 
be large. (Note that by construction, |A| < 1, and it is defined 
for a given region). Conversely, cold regions will have large, 
negative A(V). In panel (a) of Fig. 1, we have plotted A for 
regions of size 5000 km s" 1 along a mock spectrum of a simula- 
tion 7 . The mock spectrum has a jump in 7b from 7b = 2.2 x 10 4 
K below redshift z = 3.3 to 7b = 1.5 x 10 4 K above it. As ex- 
pected, there is a corresponding jump in A at z = 3.3, from 
positive values at low z to negative ones at higher z. 

To determine whether a high value of | A| is statistically sig- 
nificant - and hence whether we have identified a region with 
an unusual temperature - we proceed as follows. We repeat 
the procedure with spectra where the absorption lines are ran- 
domly scrambled (with replacement, see Theuns et al. 2002), 
in order to wash-out any correlations in the wavelet amplitudes, 
resulting from intrinsic temperature fluctuations. Using these 
randomized spectra, we can construct the statistical probability 
P(A) from the fraction of regions in the randomized spectra, 
that have a given value of A. Given P(A), we can determine 
how likely a value of A - and hence of a temperature devia- 
tion - in the original spectrum is. Performing this analysis, we 
find that the detected change in A in the mock spectrum has a 
statistical significance of more than 99.5 per cent, shown as the 
full line in panel (a) (positive values refer to regions hotter than 
average, and vice versa for negative values). This means that 
in only one of out 200 random realisation, do we, by chance, 
get values of |A| ^ 0.2. Note that we only use the simulation 
to generate the spectrum, and not to assign the statistical sig- 
nificance of a change in A. In the preparation of the mock 
spectrum, we have imposed the same biases as are present in 
the real data, by adding noise and instrumental broadening to 
the simulated lines, and by scaling the mean absorption to the 
observed value. We therefore believe that the method can be 
applied to real data with confidence. 

We have applied the same wavelet analysis to high-resolution 
echelle spectra of quasars Q0055-2169 (emission redshift z em = 
3.6, Kim et al 2001), the combined spectrum of Q0302-003 (z = 
[3 , 3 .27] , Kim et al 200 1 ) and APM 0827+5255 (z = [3 .27 , 3 .7] , 
Ellison et al 1999) and Q1422+231 (z em = 3.6, Rauch et al 1997) 
(panels b, c and d of Fig. 1 respectively). In each spectrum, 
we find a cold region at high redshift, and a hot region at lower 
redshift, each significant at the more than 99 percent level when 
compared to randomized spectra. Panels (b) and (c) appear very 

7 The cosmological parameters for this vacuum-energy dominated, flat, 
(Cl m ,Cl A ,Cl b h 2 ,h,a s ) = (0.3,0.7,0.019,0.65,0.9). The simulation box is 12/T 1 N 



similar to the mock spectrum of panel (a), where we had im- 
posed a sudden temperature increase below z = 3.3. Note that 
the implied temperature evolution is exactly opposite of what 
one would expect from photo-heating in the optically thin limit, 
in which case the IGM will gradually cool down. We take this 
as strong evidence that a large fraction of He II is reionized 
around redshift z ~ 3.3. 




redshift 

FIG. 1. — Temperature statistic A (right-hand scale, dot-dashed line) 
and statistical significance for A in per cent (left-hand scale, full line) as a 
function of redshift for the mock spectrum, Q0055-2169, the combined spec- 
trum of Q0302-003 (z = [3, 3.27]) and APM 0827+5255 (z = [3.27,3.7]), and 
Q1422+231, panels a-d respectively. Negative values for A and the signifi- 
cance denote cold regions, and vice versa for hot regions. The mock spectrum 
has a jump in T from T = 2.2 x 10 4 K below z = 3.3 to T = 1.5 x 10 4 K 
above z = 3.3. The statistical significance of the regions with a large value of 
A | at high and low z, when compared to randomized spectra made from this 
spectrum, but without a sudden temperature change, is more than 99 percent. 
The observed data exhibit a similarly significant jump, in the sense that the 
higher redshift halves are colder than the lower redshift halves. 



The mean wavelet amplitude scales approximately inversely 
with temperature: (A) oc r o _1 (Theuns et al. 2002). We have 
used hydrodynamical simulations to calibrate the proportional- 
ity constant, and investigate its dependence on 7. Applying the 
calibration to the above QSO sample, we obtain values for 7b 
in good agreement with those obtained by Schaye et al (2000) 
from the cutoff in the line width - column density relation, i.e., 
7b w 10 4 1 for z > 3.4, and T () w 10 4 ' 3 at z « 3.0. This assumes 
the value of 7 as determined by Schaye et al., but reasonable 
changes of 7 do not change 7b by more than ~ 10 per cent (The- 
uns et al. 2002). Simulations that include radiative transfer are 
required to investigate whether such a temperature change is 
consistent with He II reionization. 

Other heating mechanisms, for example shock heating by 
galactic winds, do not have a major influence on the value of T Q 
deduced from fitting the cutoff in the b-Nm diagram, at least 
as long as the volume fraction of shocked gas remains small. 
This is because the method is based on identifying the narrow- 
est lines in a region, irrespective of whether there is also a set of 
much broader lines. In contrast, the wavelet method used here 
examines all lines in a stretch of spectrum. So the fact that both 
methods find similar values for Tq suggests that photo-heating 
is indeed the dominant heating mechanism, and that the volume 

cold dark matter, smoothed particle hydrodynamics (SPH) simulation are 
□c on a side, and gas and dark matter are represented with 256 3 particles each. 



Theuns et al. 
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FIG. 2. — Temperature evolution of the IGM after z = 3.4. The points with 
error bars are determined from fitting the cutoff in the 6-Sm relation (Schaye 
et al 2000). The curves emanating from the star all assume that Tq = 10 4 4 K 
at z = 3.4 as a result of Hell reionization, and differ in the assumed UV- 
background at lower redshifts. The solid line includes photo-heating by H I, 
He I and He II. The short dashed curve neglects heating from He II, the long 
dashed curve neglects heating from both H I and He I, and finally the dot- 
dashed curve ignores photo-heating all together. All models fit the inferred 
temperature evolution extremely well. 

filling factor of gas that has been shocked by winds is small. 

2.2. Thermal evolution 
After reionization, the evolution of Tq is given by 
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where H is the Hubble parameter, kg Boltzmann's constant, 
[i the mean molecular weight, and A e is the effective radia- 
tive cooling rate [in ergg _1 s -1 ]. A e is negative (positive) for 
net cooling (heating) and includes photo-electric heating and 
cooling via recombination, excitation, inverse Compton scatter- 
ing, collisional ionization, and bremsstrahlung (we use the rates 
listed in Appendix B of Theuns et al. 1998). For gas around the 
mean density, the dominant cooling process is the adiabatic ex- 
pansion of the universe (the first term on the right-hand side of 
equation [j]), except at z > 7 where inverse Compton cooling off 
the cosmic microwave background is more efficient. The radia- 
tive heating and cooling rates depend on the ionization balance 
of the gas, which itself depends on the temperature. By cou- 
pling equation [l] to the differential equations for the ionization 
balance (listed in Appendix B of Theuns et al. 1998), we can 
thus solve for the evolution of To given a model for the ioniz- 
ing background, an initial value for Tq, and the initial ionization 
state. We have tested this procedure using full hydrodynamic 
simulations and find that the evolution of To is reproduced very 
well. 

We assume the universe to be permeated by a uniform UV- 
background with a power-law spectral shape, 
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where is the H I ionization threshold. (Note that we al- 
low for a step between the intensities of H I and He II ionizing 
photons, but normalize Jneii at the hydrogen ionization edge.) 
We set /hi = 4 x 10~ 23 , which yields a photo-ionization rate of 
Thi = 10~ 13 s -1 and vary Jueii (and the corresponding ioniza- 
tion rate Teen) as described below. Because the photo-heating 
rate is independent of the amplitude of the UV-background as 
long as the gas is highly ionized, the exact value of Jm after 
hydrogen reionization is unimportant. Since He II is not always 
highly ionized, the thermal evolution does depend on /neii- 

The predicted evolution of Tq is compared with the data from 
Schaye et al (2000) in Figure 2. We assume the IGM to be 
in ionization equilibrium at temperature Tq = 10 4 4 K at red- 
shift z = 3.4 as a result of He II reionization. The curves em- 
anating from the big star then show the subsequent evolution 
of To, for the following imposed UV-backgrounds. The solid 
curve is for a constant ionizing rate of Thi = THeii = 10~ 13 s , 
and the corresponding photo-heating from all species. Increas- 
ing or decreasing Thi by a factor 10 does not change the evo- 
lution appreciably. The short and long dashed curves ignore 
photo-electric heating of He II and of both H I and He I re- 
spectively. Finally, the dot-dashed line ignores photo-heating 
all together. The measurements of Schaye et al. (2000) below 
z ~ 3 are clearly consistent with a thermal evolution dominated 
by adiabatic cooling. 

Prior to He II reionization, some regions will already be ion- 
ized in He II by local sources. If most of the universe is reion- 
ized significantly later, then such differences in reionization 
epoch will lead to spatial variations in To- The wavelet anal- 
ysis by Theuns et al. (2002) can detect variations in To of order 
50 per cent, over a region of size 5000 km s" 1 , yet no such 
fluctuations were found in the data. The wavelet analysis by 
Zaldarriage (2001) also failed to detect any such temperature 
fluctuations. This suggests that the temperature increase is the 
result of the overall increase in the far UV-background follow- 
ing the percolation of He III regions, which prior to reionization 
are too small to be detected by current methods. 



3. HYDROGEN REIONIZATION 

Given the success in reproducing the z < 3.4 temperature 
evolution, we now turn to higher redshifts. Abel & Haehnelt 
(1999) performed radiative transfer calculations of the expan- 
sion of an ionization bubble around a QSO and found typical 
post-reionization temperatures of 4 x 10 4 K, including He II 
reionization. To put a conservative upper limit on the hy- 
drogen reionization redshift zm, we assume a higher value of 
Tq = 6 x 10 4 K. Assuming the gas to be in ionization equilib- 
rium just after reionization, we compute 7b(z) for a given zm 
as before. In Fig. 3 the predicted thermal evolution for various 
values of zhi is compared with the data. 

The temperature decreases rapidly with decreasing redshift 
because at these temperatures photo-heating cannot compen- 
sate for the rapid expansion cooling. Consequently, the tem- 
perature drops quickly below the values measured at z ~ 4, 
unless hydrogen reionization occurred relatively recently. The 
solid curves are for Thi = 10~ 13 s" 1 and a range of reioniza- 
tion redshifts, as indicated in the figure. For zh = 6.2, we show 
Thi = 10~ 14 s -1 for comparison. The models with late H I reion- 
ization, zh < 8, fit the data best. Models with Zh > 9 have z ~ 4 
temperatures that are significantly below the measured values. 
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For example, the zh = 9.2 curve has a reduced \ 2 of 4.5 for the 
5 z > 3.4 data points. 
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FIG. 3. — The temperature evolution of the IGM above redshift 3.4. The 
solid curves indicate the evolution of the temperature at the mean density for 
various HI reionization redshifts zh, as indicated. The post-hydrogen reion- 
ization temperature is assumed to be Tq = 6 X 10 4 K and the hydrogen photo- 
ionization rate is Thi = 10~ 13 s (the short dashed line has Thi = 10~ 14 s ). 
The He II photo-ionization rate is adjusted so that the He III abundance 
*HeM ~ 0.1 at z = 3.5. The solid line connecting filled squares is for zh = 10.2, 
and a higher Hell photo-ionization rate, .th c iii(z = 3.5) = 0.6. Finally, the 
long dashed line has zh = 20, but a still higher He II photo-ionization rate, 
*Hein(z = 3.5) = 0.95. If He is mostly singly ionized at z > 3.5, then the rapid 
decrease in Tq after reionization places an upper limit of z.h < 9 on the redshift 
of hydrogen reionization. 

The contribution of He II photo-heating at z > 5 is uncer- 
tain. It is likely that stars dominate the H I ionizing back- 
ground at such high redshifts (e.g., Madau, Haardt, & Rees 
1999), and stars emit very few He II ionizing photons. How- 
ever, very massive and extremely metal poor stars could lead to 
a non-negligible He II ionizing background (Tumlinson & Shull 



2000). The QSO contribution is also uncertain, although the 
paucity of faint point sources in the Hubble Deep Field does 
provide some constraints (Haiman, Madau & Loeb 1999). If 
He II reionized at z ~ 3, as argued in the previous sections, then 
the He III abundance, XHein, should be small at higher z- We 
therefore set Teen = 10~ 15 s" 1 so that XH e m ~ 0.1 at redshift 
z = 3.5. This limits the plausible contribution from He II photo- 
heating, and allows us to put a conservative upper limit on the 
redshift of H I reionization, zh < 9. 

The importance of He II heating is illustrated by the solid 
curve connecting the filled squares. This model has zh = 10.2, 
but the He II ionizing background is increased such that JtHeiii ~ 
0.6 at z= 3.5. This increases Tq(z = 4) significantly, although it 
still falls below the measured values. Finally, the long dashed 
line is for a model with zh = 20 and Jueii further increased so 
that XHein ~ 0.96 at z = 3.5. The temperature of this model is 
consistent with the data, yet He II is ionized at the more than 90 
percent level as early as z = 5. Such a high level of ionization 
conflicts with the evidence that He II reionizes at z ~ 3 from the 
observed He II opacities and the associated increase in 7q. 

All this leads us to the following conclusions. Two inde- 
pendent methods consistently find a rather sudden increase in 
the temperature of the IGM over the range z ~ 3.5-3.0, which 
we associate with He II reionization. If this interpretation is 
correct, then the He III fraction must be low at higher redshifts. 
Therefore, above redshifts 3.5, He II photo-heating cannot be 
significant and the IGM cools rapidly following H I reioniza- 
tion. The high values of the IGM temperature at z ~ 4 then 
require that H I reionization occurred late as well, zh < 9, 
for any reasonable value 7b <J 6 x 10 4 K for the HI reioniza- 
tion temperature. More plausible reionization temperatures of 
7q ~ 4 x 10 4 K and 2 x 10 4 K would constrain the hydrogen 
reionization redshift further to zh < 8 and zh < 7, respectively. 
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